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I.  INTRODUCTION 


r 


Combustion  processes  normally  involve  many  chemical  species,  related  by 
complex  kinetics  and  often  strongly  influenced  by  molecular  transport  and 
fluid  flow.  The  difficult  problem  of  combining  the  chemistry  and  transport  is 
commonly  made  tractable  by  idealizing  the  complex  kinetic  network  as  a  single 
overall  or  "global"  reaction.  Such  a  strategy  is  potentially  useful  when  the 
detailed  chemical  steps  are  unknown  or  when  the  difficulty  of  the  flow  problem 
necessitates  simplified  chemistry.  Relatively  little  work  has  been  done, 
however,  to  assess  the  fidelity  with  which  a  global  reaction  can  represent  the 
net  effects  of  the  complex  reaction  network.  This  question  is  addressed  here 
for  the  particular  case  of  premixed,  laminar,  one-dimensional,  steady-state 
flames. 


Previous  work  in  this  vein  was  reported  by  Levy  and  Weinberg,  who 
obtained  temperature  profiles  through  very  lean  (and  very  slow)  ethylene-air 
flames  by  the  ray  deflection  method.  Using  a  form  of  the  energy  conservation 
equation,  they  transformed  the  temperature  profile  into  a  profile  of  the 
volumetric  rate  of  heat  release  through  the  flame  front.  The  authors 
associated  this  net  rate  of  heat  release  with  that  resulting  from  a 
hypothetical  global  reaction.  A  flame  speed  was  then  calculated  from  the 
"measured"  reaction  rate  profile  by  means  of  the  single-reaction  flame  theory 
of  Spalding. ^  Close  agreement  was  noted  between  this  calculated  result  and 
the  directly  measured  flame  speed  for  a  (rather  narrow)  range  of 
stoichiometries  and  initial  temperatures.  Such  a  finding  gave  encouragement 
to  the  notion  that  the  flame  propagation  formalism  based  on  a  single  reaction 
might  indeed  have  some  utility. 


In  a  subsequent  paper  Levy  and  Weinberg-^  probed  the  overall  kinetics 
concept  directly.  By  examining  the  heat  release  rates  at  the  same  temperature 
for  a  number  of  different  flames,  they  determined  that  the  effective  reaction 
orders  (with  respect  to  fuel  and  oxidizer  concentrations)  vary  greatly  through 
the  flame.  Using  the  heat  release  profiles,  they  also  found  that  the 

>  effective  activation  energy  typically  increases  with  temperature  through  the 

>  flame  front  then  decreases  near  the  maximum  temperature.  The  conclusion  was 

j  that  no  single  reaction  (of  an  Arrhenius  form)  could  account  for  the  thermal 

1  structure  of  the  flame.  However,  they  did  note  that  the  values  of  apparent 

j  activation  energy  and  total  order  near  the  maximum  heat  release  rate  were 


1.4 .  Ijp.vy  and  E.J.  Weinberg,  "Optical  Flame  Structure  Studies:  Some 
Conclusions  Concerning  the  Propagation  of  Flat  Flames,"  7th  Intern.  Symp. 
Combustion,  pp.  296-303,  1959. 

O 

D.B.  Spalding,  "I.  Predicting  the  Laminar  Flame  Speed  in  Gases  with 
Temperature-Explicit  Reaction  Rates ,"  Combustion  and  Flame,  Vol.  1, 
pp.  287-295,  1957.  "II.  One-Dimensional  Laminar  Flame  Theory  for 
Temperature-Explicit  Reaction  Rates ,"  Combustion  and  Flame,  Vol.  1, 
pp.  296-307,  195  7. 

7  t 

A.  Levy  and  E.J.  Weinberg,  " Optical  Flame  Structure  Studies:  Examination  of 
Reaction  Rate  Laws  in  Lean  Ethylene-Air  Flames,"  Combustion  and  Flame, 

Vol.  3,  pp.  229-253,  1959. 
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close  to  those  obtained  at  blowout  in  a  stirred  reactor  experiment.  This  last 
observation  seemed  to  justify  use  of  the  overall  reaction  concept  for  those 
phenomena  strongly  dependent  on  the  maximum  energy  release  rate. 


Our  objective  in  this  paper  is  similar  to  that  of  Levy  and  Weinberg; 
however,  here  we  make  use  of  a  laminar  flame  model  which  includes  multispecies 
ami  mu  1 1 i reaction  chemical  kinetics.  This  approach  offers  the  advantage  of 
complete  information  about  the  system  as  well  as  access  to  a  wide  range  of 
flames  and  stoichiometries.  The  only  potential  drawback  lies  in  our  uncertain 
knowledge  of  the  many  reaction  steps  and  rates  for  some  flames.  Because  our 
aim  is  heuristic,  this  trait  is  not  considered  a  serious  limitation.  Our 
procedure  is  to  solve  the  flame  equations  with  the  complex  kinetics,  obtain  an 
overall  rate  for  the  flame  from  the  heat  release  profile,  solve  the  flame 
equations  for  this  single  reaction,  and  then  make  detailed  comparisons  between 
the  two  solutions.  This  is  done  for  a  wide  range  of  stoichiometries  of 
several  flames. 

Our  results  suggest  that  a  single  reaction  formalism  does  a  credible  job 
ot  representing  the  thermal  structure  (and  even  the  reactant  profiles)  of 
complex  flames.  This  conclusion  is  not  as  contradictory  to  the  Levy  and 
Weinberg  study^  as  it  may  seem.  Their  conclusion  was  based  on  a  rather 
restrictive  diagnosis  of  the  heat  release  profile.  A  comparable  analysis  of 
our  results  leads  to  similar  findings,  i.e.,  the  activation  energy  varies 
through  the  flame  front.  Our  conclusion  follows  from  the  relative 
insensitivity  of  thermal  structure  to  the  reaction  rate  details.  This  can 
only  be  seen  by  actually  solving  the  flame  equations  using  an  appropriate 
overall  reaction. 

Westbrook  and  Dryer^  have  developed  a  procedure  for  determining  global 
kinetics  parameters  from  the  experimental  flammability  limits  and 
stoichiometric  flame  speed  at  atmospheric  conditions.  A  single-reaction  flame 
model  using  these  parameters  can  be  expected  to  predict  reasonable  flame 
speeds  over  a  wide  range  of  equivalence  ratios  from  very  lean  to  very  rich. 

The  procedure  may  prove  useful  to  those  applications  in  which  flame  speed  is 
of  decisive  importance.  However,  an  application  of  interest  to  us  is 
propellant  combustion  where  two-phase  flame  holder  effects  give  importance  to 
temperature  gradients  near  the  flame  front.  Calibration  of  the  global 
parameters  by  flame  speeds  alone  does  not  guarantee  accurate  temperature 
profiles,  but  both  the  flame  speed  and  temperature  profile  can  be  computed 
from  an  accurate  heat  release  profile.  The  thrust  of  this  paper  is  that  one 
can  find  a  global  reaction  which  mimics  the  original  heat  release  profile 
sufficiently  well  that  accurate  flame  speeds  and  temperature  profiles  (at 
least  through  the  flame  front)  can  be  computed  from  the  global  description. 

The  significance  of  this  conclusion  is  not  necessarily  diminished  by  the 
finding  that  the  optimum  global  kinetics  values  vary  with  stoichiometry.  In 


4  ^ 

<  • K .  Westbrook  and  F.L.  Dryer,  "Simplified  Reaction  Mechanisms  for  the 
Oxidation  of  Hydrocarbon  Fuels  in  Flames,"  Comb.  Sci.  Tech.,  Mol.  27, 
pp.  31-43,  1981.  - 
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application  work,  generality  might  be  achieved  simply  by  parameterizing  the 
kinetics  values  in  terms  of  stoichiometry,  pressure,  etc.,  as  suggested  by 
Dryer  and  Westbrook. 

Our  interest  in  the  global  reaction  concept  springs  from  two  potential 
application  schemes.  The  first  assumes  that  the  detailed  multireaction 
kinetics  are  known,  but  so  complex  as  to  preclude  their  use  in  practical  flow 
problems.  One  might  then  determine  the  appropriate  global  parameters  in  a 
manner  similar  to  that  used  in  this  paper.  The  alternative  scheme  assumes 
that  little  is  known  of  £he  detailed  kinetics  of  a  flame.  It  is  for  this  case 
that  Westbrook  and  Dryer"4  developed  their  method  of  predicting  flame  speeds 
vs.  equivalence  ratio  from  only  a  few  selected  measurements.  To  realize  the 
thermal  structure  as  well  as  the  flame  speed,  more  extensive  measurements  are 
required  to  calibrate  the  global  reaction.  The  procedure  would  involve 
measuring  temperature  profiles  in  the  subject  flame  using  fine  thermocouples 
or  laser  Raman  techniques  and  then  computing  heat  release  profiles  via  the 
energy  conservation  equation.  ’  ’  Determination  of  the  global  parameters 
would  then  proceed  as  described  in  this  paper.  Prerequisite  to  either  of 
these  two  approaches  is  the  assurance  that  global  reactions  can  reproduce  the 
thermal  structure  of  real  flames  with  reasonable  fidelity.  This  paper 
provides  both  this  assurance  and  a  systematic  procedure  for  determining 
appropriate  global  reaction  parameters  from  a  heat  release  profile,  which  is 
either  calculated  from  known  detailed  kinetics  or  deduced  from  a  measured 
temperature  profile. 


II.  DETAILED  CHEMISTRY  MODEL 

In  a  previous  paper,®  the  equations  describing  a  one-dimensional, 
laminar,  premixed,  steady-state  flame  were  derived  for  different 
approximations  to  the  multicomponent  transport  coefficients.  We  will  use  tin 
simplest  approximat ion  (Method  V),  which  still  gives  accurate  results.  The 
corresponding  equations  are  the  following. 

Conservation  of  energy  is  expressed  by 


c 

P 


m 

o 


dT  d_ 
dx  +  dx 


(\  ~  )  -  E  R.  M. 
dx  11 


h.  =  0 


( n 


V.L.  Dryer  and  C.K.  Westbrook,  NATO  AGARD  Conference  Proceedings  No. 
pp.  1 4-1  to  14-1  7,  1979. 

ft 

R.  Pruhimm  < md  E.  Rurke,  "Measurement  of  Temperature  Distribution  in  a  Loi.'- 
Pressure  Flat  Flame,"  J.  Chem.  Phys.,  Vol.  22,  pp.  824-830,  1984. 

'J.A.  Hicks,  "The  Low-Pressure  Decomposition  Flame  of  Ethyl  Nitrate 8th 
Symposium  (International)  on  Combustion,  The  Combustion  Institute,  1982, 
pp.  487-496. 

q 

T.P.  Coffee  'md  9.M.  Heimerl,  "Transport  Algorithms  for  Premixed,  leminar 
Steady-State  Flames,"  Combustion  ani  Flame,  Vol.  43,  pp.  273-289,  1981. 
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where  T  is  the  temperature  of  the  mixture,  c  is  the  specific  heat  of  the 
mixture,  m0  is  the  mass  flux  through  the  flame,  \  is  the  thermal  conductivity 
of  the  mixture,  R^  is  the  rate  at  which  species  I  is  produced  (or  consumed)  by 
chemistry,  Mj  is  the  molecular  weight  of  species  i  and  hj  is  the  specific 
enthalpy  of  species  i. 

The  species  continuity  equations  are  given  by 


dY. 

1 

n  j — 
o  dx 


,  dY 

2-  (p  D.  ~r  )  +  R.  M  =  0, 
dx  tm  dx  i  i 


(2) 


where  Y^  is  the  mass  fraction  of  species  i,  P  is  the  density  of  the  mixture, 
and  r»lm  is  the  diffusion  coefficient  of  species  i  into  the  mixture.  The 
boundary  conditions  are 


'  and  Y. 

=  Y  , 

i  =  1,  2,  .. 

u  i 

iu 

dY, 

4  dT 
and  -r—  = 
dx 

_ i  _ 

dx 

0,  i  =  1,2,. 

at  x  = 


In  the  above  equations,  the  effects  of  radiation,  viscosity,  and  body 
forces  are  ignored.  Since  the  burning  velocity  is  smalL  compared  with  the 
local  speed  of  sound,  the  pressure  is  taken  to  be  constant.  Besides  these 
standard  assumptions,  we  have  also  assumed  (see  Ref.  8  for  discussion)  that 
thermal  diffusion  is  negligibly  small, ^that  Fick's  law  holds  for  the  diffusion 
velocities,  and  that  the  quantities  p  ,  P^ ,  and  c  are  constant.  The 
constants  are  chosen  a  priori.  m  ^ 


The  equations  above  are  for  the  steady  state  solution.  The  numerical 

procedure  actually  uses  the  time  dependent  equations  and  a  transformed  space 

coordinate.  Starting  with  arbitrary  profiles,  the  equations  are  integrated  in 

time  until  the  steady  state  solution  is  obtained.  The  burning  velocity, 

S,  of  the  flame  is  defined  as  the  velocity  of  the  flame  relative  to  the  fluid 

at  rest,  i.e.,  at  x  =  -  °°.  Since  in  our  coordinate  system  the  flame  is  at 

rest,  S  =  v(-°°)  =  m  /p  .where  p  is  the  density  of  the  unburned  mixture. 

o  u  u 

The  production  rate  functions,  R^(T),  are  defined  by 
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\,r>  Ikfr 
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T.P .  Coffee  and  J.M.  Heimerl,  "A  Method  for  Computing  the  Flame  Speed  for  a 
Laminar ,  Premixed,  One  Dimensional  Flame,"  BRL  Technical  Report,  ARBRL-TR- 
02212,  January  1980  (AD  A082803). 

1°T.P.  Coffee,  "A  Computer  Code  for  the  Solution  of  the  Equations  Governing  a 
Laminar,  Premixed,  One-Dimensional  Flame,"  BRL  Memorandm  Report,  A RB RL -MR- 
03 165,  April  1982  (ADA114041). 
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where  Vp  r  and  vi  r  are  the  number  of  molecules  of  species  i  entering  into 
reaction  r  as  product  and  reactant,  respectively.  NR  is  the  number  of 
reactions,  Cj  is  the  concentration  of  species  j  (pY./M.),  and  kpr  and  k^r  are 
the  forward  and  reverse  rate  coefficients,  respectively!  The  rate 
coefficients  depend  on  the  temperature  through  the  generalized  Arrhenius 
relat  ions 

B 

k_  =  A,  T  exp(-  Ef  /RT)  (6) 

tr  fr  fr 

and 

khr  '  V  r‘'’re'<’  <  -  V/RTK  <7> 

where  K  is  the  gas  constant,  E  is  the  activation  energy  and  A  is  the  frequency 
factor  for  the  reaction. 

The  net  rate  of  heat  release  per  unit  volume  is  defined  by 

q  =  -  £  R.M.h.,  (8) 

111 

It  is  through  this  term  ( Eq .  1)  that  the  chemical  kinetics  influences  the 
temperature  profile. 


III.  ONE  REACTION  MODEL 

We  will  now  assume  a  single  irreversible  overall  reaction  of  the  form 
fuel  plus  oxidizer  goes  to  product,  that  is, 

v  F  +  v  0  +  v  p,  (9) 

fop 

where  Up,  vQ,  and  v  are  the  stoichiometric  coefficients.  The  governing  Eqs. 
(1)  through  (4)  remain  the  same.  Equation  (5)  can  be  simplified.  We  obtain 


RF  -  -  "F  k  (CF>  F  <c0>  °> 

Ro  ’  '  k  <V  F  <Co> 


and 


\)  v 

R  =  y  k  (CJ  F  (c  )  ° 

P  p  F  O 


(  10) 

(11) 

(12) 


The  heat  release  is  then 


<1  =  MfV  (CF)  (Co) 


>  M  v 

O  .  o  o 

Ihf  V7 


M  v 

h  - 

o  Mpv  p 


(13) 


1  ) 


We  define 


m  v 
no  o 


^  hF  +  M.v 


F  F 


M  v 
P  P 

M  v 
F  F 


(14) 


whore  Q  is  the  heat  of  reaction  per  unit  mass  of  fuel. 

O 

From  the  therinodynami cs  the  relation 


c  (T.  -  T  )  =  l  h  (Y.  -  Y..  ) 

p  b  u  i  in  lh 


(15) 


can  he  obtained,  where  is  the  temperature  of  the  completely  burned  mixture 
and  the  Yiv  are  the  mass  fractions.  Since  there  is  only  a  single  reaction, 
however,  any  change  in  the  fuel  must  be  matched  by  a  change  in  the  oxidizer 
and  product;  that  is. 


M  v 

y  -  Y  =  — — —  ( Y  -  Y  ) 

on  oh  Mv  Fu  Fh 


and 


Y 

pu 


ph 


M  v 

-fi-E  (Y 

Vf  pu 


Y  ,). 
pb 


(lb) 


(17) 


Substituting  into  (15)  and  solving  for  Q, 


Q  =  c 


(Tb" 


T  )  / 
u 


(Y 


Fu 


V 


(IK) 


The  reaction  rate  k  is  expressed  as  an  unmodified  Arrhenius  rate 

k  =  exp  (-E/RT) .  ( 19) 

For  convenience,  the  constants  in  Eq .  (13)  will  be  absorbed  into  a  new 
constant  A,  giving 

v  v 

q  =  Q  (P  Y„)  F  (PY  )  °  A  exp  (-E/RT).  (20) 

r  O 

Equation  (20)  indicates  that,  given  a  solution  of  the  detailed  chemistry 
model  (q),  a  set  of  constants  Q,  A,  and  E  may  be  found  for  which  the  one 
reaction  model  most  closely  follows  the  detailed  model. 


IV.  DERIVING  THE  ONE  REACTION  MODEL  PARAMETERS 

In  this  section  we  derive  expressions  for  the  one  reaction  model  such 
that  Eq.  (20)  can  be  used  for  finding  the  kinetic  parameters,  given  a  heat 
release  profile  q  from  the  detailed  chemistry  model. 
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Equation  (20)  expresses  q  as  a  function  of  T,  Yp,  YQ,  and  Yp,  since 


However,  it  is  inconsistent  to  use  the  values  of  Yp,  YQ,  and  Y  from  the 

detailed  chemistry  model.  In  the  detailed  model,  there  are  otner  species,  so 

that  Y„  +  Y  +  Y  *  1 .  In  the  simplified  model,  by  assumption  there  are  no 
F  .  o  p 
other  species. 


Consequently,  we  make  an  additional  assumption  that  all  the  Lewis  number 
are  one;  that  is. 


P  D .  c 
im  p 

X 


(22) 


for  fuel,  oxidizer 
good  approximation 


d_ 

dx 


and  product. 


Eqs. 

(i) 

and 

m 

i  o 

dT 

dx 

-  Q 

For  all  except  the  light  species,  this  is 
(2)  can  he  written 


RpHp  =  -  q 


(23) 


a 


and 


dYF  ,  dYF 

- —  )  -  m  -t — 
dx  o  dx 


RpMp  . 


(24) 


The  equation  in  terms  of  (-cp  T/Q)  and  Yp  are  formally  identical.  Also 
matching  boundary  conditions,  we  obtain 


\  ~  YFb  *  %  "b  -  T>/Q-  <25) 

So  given  T,  Yp  can  be  computed,  and  also  YQ  and  Y  since  the  stoichiometric 
coefficients  are  known  for  the  overall  reaction,  Eq.  (9). 

The  heat  release  can  now  be  viewed  as  a  function  only  of  T.  Moreover, 
the  one  reaction  model  simplifies  to  the  single  Eq.  (23). 


Tlu*  heat  of  reaction  i)  can  he  computed  from  Kq .  (14),  since  the 
enthalpies  are  known  quantities.**’*-  i)  may  also  he  computed  from  Kq  .  (18). 
The  adiabatic  or  burned  temperature  and  the  mixture  specific  heat  r  can  be 
computed  from  thermodynamic  considerations.  Since  the  one  reaction  i n  the 
simplified  model  is  assumed  irreversible,  Yj..^  can  be  computed  knowing  Yj.-U  and 
YiU  .  So  all  the  quantities  on  the  right  side  of  Kq .  (18)  are  known. 

However,  Kqs.  (.14)  and  (18)  are  Inconsistent.  Thi  s  occurs  because  the 
two  models  have  different  adiabatic  species  concentrations.  In  particular, 
some  intermediates  ami  radicals  wilt  still  exist  at  the  burned  conditions. 
And,  since  the  reactions  are  reversible,  combustion  may  ho  incomplete.  In 
order  to  preserve  the  proper  adiabatic  temperature,  Kq .  (18)  will  he  used  to 
compute  i). 


What  renut  l ns  is  to  find  "best  fit"  values  ol  A  and  K  in 


v  v 

q/i)|  (pY  )  (pY  )  °|-  A  exp  (-K/KT)  . 

r  O 


(2b) 


Tlie re  can  he  a  difficulty  here.  The  expression  on  the  right  will  ho 
monotoni ca 1 lv  increasing  for  positive  A  and  K.  Kor  the  cases  considered,  the 
expression  on  the  left  normally  reaches  a  peak  somewhere  past  the  q  peak  and 
then  decreases.  We  then  only  fit  the  expression  up  to  the  peak.  This  means 
the  post  flame  region  will  not  be  well  resolved.  However,  tills  region 
probably  will  not  he  well  resolved  by  any  one  reaction  model.  The  important 
reactions  in  this  region  are  radical  recombinations,  and  radicals  are  ignored 
by  the  simplified  model.  The  fit  is  m^<jc  using  an  absolute  error  criterion 
with  proper  transformation  of  weights. 

In  using  the  fitted  values  of  A  and  K  to  solve  Kq .  (2)),  we  must  be 
concerned  with  the  cold  boundary  difficulty.  That  is,  from  Kq .  (20),  heat  is 
generated  at  any  temperature,  including  the  temperature  of  the  unburned 
mixture  T  .  (liven  enough  time,  the  mixture  will  spontaneously  combust.  This 
contradicts  our  idealization  of  an  infinite  unreacting  column  of  gas  with  the 
flame  front  propagating  through  the  mixture.  In  practice,  we  will  numerical  lv 
integrate  over  a  finite  interval.  As  long  as  the  heat  release  at  T  is  small, 
the  heat  generated  before  the  gas  encounters  the  flame  front  will  he 
negligible.  This  is  true  for  most  of  the  cases  studied  in  this  paper. 

However,  for  a  few  cases  the  heat  release  is  too  large  at  T  ,  and  the  cold 
boundary  conditions  will  he  distorted. 
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To  avoid  this,  we  utilize  the  concept  of  an  ignition  temperature.*^1  For 
T  <  T  +  10  K,  the  heat  release  q  is  set  to  zero.  For  T  >  Tu  +  10  K,  F.q.  (20) 
is  used.  That  is,  the  gas  cannot  ignite  until  it  reaches  the  ignition 
temperature  Tu  +  10.  This  circumvents  the  cold  boundary  difficulty,  while 
having  little  effect  on  the  rest  of  the  flame. 

So  far  we  have  assumed  that  v  and  vq  are  known.  It  is  also  possible  to 
make  a  3  parameter  fit,  where  v p  is  a  new  parameter  to  be  chosen  and  the 
relation  between  Vp  and  v  is  determined  by  the  overall  stoichiometry.  Then 
the  quantity  vp  +  v  is  the  overall  order  of  the  reaction.  We  will  also 
discuss  this  type  of  fit. 

There  are  several  differences  between  the  above  procedure  and  that 
recommended  by  Westbrook  and  Dryer. 4  The  primary  distinction  is  that 
Westbrook  and  Dryer  fit  flames  speeds,  while  we  fit  the  heat  release. 

Correctly  fitting  the  thermal  structure  will  lead  to  the  proper  flame  speed; 
the  converse  is  not  true.  Also,  the  Westbrook  and  Dryer  model  leads  to  an 
incorrect  adiabatic  temperature;  since  we  are  interested  in  the  thermal 
structure,  our  procedure  has  been  developed  so  as  to  reproduce  the  proper 
adiabatic  temperature.  Since  we  assume  unit  Lewis  numbers,  our  global  model 
reduces  to  only  t lie  temperature  equation,  while  the  Westbrook  and  Dryer  model 
includes  equations  for  the  reactants  and  products.  Also,  in  the  Westbook  and 
Dryer  model,  a  single  set  of  parameters  is  selected  which  fit  a  range  of 
stoichiometries.  As  we  shall  see,  our  parameters  will  vary  with 
stoichiometry,  reflecting  the  changes  in  thermal  structure  due  to  different 
initial  conditions. 


V.  FLAMES  TO  BE  CONSIDERED 

To  increase  the  generality  of  our  conclusions,  we  will  consider  three 
different  types  of  flames. 

The  simplest  flame  mechanism  is  that  of  an  ozone  decomposition  flame.*'’ 
There  are  three  reactions,  only  two  of  which  are  important  in  the  pre-flame 
and  flame  region: 

0_  +  M  ♦  0  +  0  +  M 

3  «-  2 

0  +  °3  t  °2  +  °2  * 

In  the  post  flame  region  the  radical  recombination  reaction 

o  +  o+  m>o2  +  o2. 

becomes  important. 


*4F.A.  Williams ,  Combustion  Theory,  Addison-Wesley,  Reading,  MA,  1965. 

.M.  Heimevl  and  T.P.  Coffee,  "The  Detailed  Modeling  of  Premixed,  Laminar 
Steady-State  Flames.  I.  Ozone,"  Combustion  and  Flame,  Vol.  39,  pp.  301-315. 
1980. 
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A  more  complicated  case  is  hydrogen/air  flames.  The  mechanism  used  is 
due  to  Dixon-Lewis , 16  and  consists  of  8  species  and  17  reversible  reactions. 

Last,  we  consider  methane/air  flames.  The  mechanism  used  is  obtained 
from  Pixon-Lewis , 1^  and  consists  of  14  species  and  34  reactions.  The 
hydrogen/air  scheme  Is  a  subset  of  this  mechanism.  While  the  mechanism  is 
intended  for  lean  and  slightly  rich  flames,  we  will  also  use  it  for  rich 
flames  in  order  to  obtain  our  comparisons. 

All  flames  are  computed  for  atmospheric  pressure.  For  the  ozone  flames, 
the  initial  temperature  of  the  gas  is  300K.  For  the  other  flames,  the  initial 
temperature  is  298K.  The  overall  reactions  for  these  mechanisms  are  taken  to 
be 


I  H.;  +  02  *  2  H20, 

and  CH,  +  2  0.,  *  ClL  +  2  HO. 

*4  Z  2  Z 

First,  we  consider  the  feasibility  of  the  overall  reaction  fit.  To  do 
this,  consider  taking  the  logarithm  of  Eq.  (26);  that  is, 

v  v 

In  | q/Q  (eYp)  F  (PYQ)  °  ]  =  *n(A)  -  K/RT,  (27) 

and  graphing  both  sides  as  a  function  of  1/T.  The  right  side  is  the  equation 
of  a  straight  line.  The  shape  of  the  left  side,  where  q  is  the  heat  release 
for  the  detailed  chemistry  model,  gives  a  quick  impression  of  the  adequacy  ot 
a  single  Arrhenius  type  reaction. 

Figure  l  gives  these  graphs  for  a  50%  Ozone  flame  (oxygen  diluent),  a 
rich  H.,/air  flame  (40%  H2),  and  a  stoichiometric  CH^/atr  flame  (9.5  %  CH^). 
The  activation  energies  of  the  least  square  tits  to  q  are,  respectively,  7.9 
kcal,  6.4  kcal,  and  29.1  kcal. 

All  three  graphs  display  a  high  activation  energy  (that  is,  a  steep 
slope)  near  the  cold  boundary.  The  Arrhenius  fits  cannot  resolve  this, 
necessitating  the  concept  of  an  ignition  temperature. 

Next,  there  is  a  long  region  where  the  slope  gradually  increases.  For 
the  0-j  and  1^  flames,  this  increase  is  very  gradual,  and  the  Arrhenius  fit  is 


1ft. 

G.  Dixon-Lewis ,  "Kinetic  Mechanism  Structure  and  Properties  of  Premixed 
Flames  in  Hydrogen-Oxygen-Nitrogen  Mixtures Proc.  /?.  Soc.  Lond.  A, 

Vol.  292,  pp.  45-99,  1979.  "  . . 

7  7.  •  , 

G.  Dixon- Lewi s,  "Aspects  oj  the  Kinetic  Modeling  of  Methane  Oxidation  in 
Flames t"  First  Specialists  Meeting  (International)  of  tfa;  Combustion 
Institute,  Vol.  1,  pp.  284-289,  1981. 
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Release  Divided  by  the  Concentrations  (bine)  Versus  the  Logarithm  of  the 
Overall  Arrhenius  Fit  (Dot) 


quite  good  over  most  of  the  temperature  range.  The  CH^  flame  exhibits  a  more 
complicated  structure,  but  there  is  still  a  fairly  large  region  near  the  peak 
heat  release  where  the  Arrhenius  fit  is  quite  good. 

Finally,  there  is  a  post-flame  region,  where  the  slope  changes  sign.  An 
Arrhenius  fit  cannot  resolve  this. 

The  results  for  the  methane/air  heat  release  profile  agree  qualitatively 
with  those  of  Levy  and  Weinberg  for  an  ethylene/air  flame.  The  activation 
energy  (or  slope)  increases  with  increasing  temperature.  Since  our  single 
reaction  parameters  are  generated  by  a  least  squares  fit  witli  an  absolute 
error  criterion,^  the  fit  is  most  accurate  in  the  region  of  maximum  heat 
release,  which  Levy  and  Weinberg  found  most  influential  to  the  burning  rate. 

It  is  of  interest  to  learn  how  accurately  these  same  parameters  reproduce  the 
thermal  structure  as  well. 

VI.  TKMPKRATltKK  AN!)  SPECIKS  PRO  FILLS 

In  this  section  we  will  compare  the  profiles  generated  by  the  detailed 
chemistry  model  with  those  generated  by  the  single  reaction  model.  The  same 
three  examples  considered  in  the  last  section  will  be  used. 

Figures  2,  3,  and  4  show  the  heat  release  and  temperature  profiles  for 
our  three  examples.  The  heat  release  profiles  from  the  overall  reaction  are 
reasonable.  However,  the  CH^/air  profile  is  especially  inaccurate  for  the 
early  part  of  the  flame,  as  expected  from  considering  Figure  1.  Nevertheless, 
all  three  temperature  profiles  are  quite  accurate  through  the  flame  front.  In 
the  post-flame  region,  the  detailed  model  shows  a  slow  temperature  rise,  due 
to  the  recombination  of  radicals.  The  overall  reaction  model,  which  ignores 
radicals,  rapidly  approaches  the  adiabatic  temperature.  Consequent ty,  the 
single  reaction  model  is  somewhat  inaccurate  in  the  post-flame  region. 

Since  we  have  assummed  unit  Lewis  number  in  the  overall  reaction  model, 
species  profiles  can  also  be  derived.  Figure  5  shows  the  comparisons  for  the 
ozone  flame.  The  agreement  is  quite  good.  The  errors  are  primarily  due  to 
the  neglect  of  the  0  profile,  rather  than  the  unit  Lewis  number  assumption. 

The  agreement  becomes  better  for  slower  flames  (less  0),  and  not  quite  as  good 
for  faster  flames  (more  0).  For  the  hydrogen/air  flame,  the  mole  fractions, 
as  shown  in  Figure  b,  are  not  in  good  agreement.  In  the  detailed  model,  there 
are  large  radical  concentrations  in  the  flame  front.  These  lead  to  rapid 
consumption  of  the  H9  and  through  the  chain  branching  reactions.  The  one 
reaction  model,  whicfi  has  a  low  activation  energy,  cannot  match  this.  For  a 
slower  flame  (either  richer  or  leaner),  the  agreement  is  much  better. 

Figure  7  shows  the  major  species  profiles  for  the  CH^/air  flame.  The 
main  errors  in  the  species  profiles  for  a  CH^/air  flame  are  due  not  to 
radicals,  which  have  small  concentrations  in  all  cases,  but  to 
intermediates.  The  overall  reaction  model  assumes  that  the  only  species  are 
CH ,  ,  C>2 ,  0(>2  >  and  H90.  There  are  several  intermediate  species,  of  which  ll9 
and  CO  are  the  most " important ,  and  except  for  lean  flames,  the  equilibrium" 
products  include  substantial  amounts  of  and  CO.  This  causes  the  009  and 
H2O  profiles  from  the  one  reaction  model  to  approach  incorrect  values."  The 
profile  agreement  is  better  for  lean  flames  and  not  quite  as  good  for  rich 
f 1 ames . 
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Figure  2.  50%  0  ,  50%  0  Flame.  Detailed  Chemistry  Model  (Line).  Overall  Reaction  Model  (Dot) 


Figure  3.  402  H  ,  602  Air  Flame.  Detailed  Chemistry  Model  (Line).  Overall  Reaction  Model  (Dot) 


Figure  5.  SOX  0  ,  50.  0  Flar.e.  Detailed  Chemistry  Mode  1  CLinej.  'Overall  Reaction 


Figure  6.  402  H  ,  602  Air  Flame.  Detailed  Chemistry  Model  (Line).  Overall  Reaction  Model  (Dot) 
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Overall,  the  species  profiles  are  reasonable,  considering  the 
simplification  in  the  one  reaction  model. 

VII.  FLAMK  Sl’KKDS 

In  this  section  we  will  consider  more  closely  the  parameters  generated  by 
the  overall  reaction  fits,  as  well  as  comparing  flame  speeds.  In  order  to 
discern  possible  patterns,  many  initial  stoichiometries  are  considered  tor 
each  of  the  three  flames. 

The  parameters  for  a  set  of  ozone  t lames  are  given  in  Table  I.  The 
•activation  energies  derived  are  between  those  of  reaction  1  (22.2  kcal)  and 
reaction  2  (4.b  kcal).  The  slower  and  cooler  the  flame  (that  is,  the  smaller 
the  initial  ()•>  concent  rat  ion) ,  the  higher  the  overall  activation  energy.  This 
can  tie  explained  by  studying  the  detailed  model.  For  a  slower  flame,  only 
small  amounts  of  0  are  produced.  The  0  that  is  created  reacts  almost 
immediately  via  reaction  2.  So  reaction  1  is  close  to  a  rate  limiting  step, 
and  the  overall  reaction  has  a  similar  activation  energy.  For  a  fast  flame, 
large  amounts  of  0  are  produced,  which  diffuses  ahead  of  the  flame.  The  flame 
propagation  is  determined  to  a  large  extent  by  the  low  activation  energy 
reaction  2.  This  is  reflected  in  the  fitted  values  of  E.  Table  2  compares 
the  flame  speeds  from  the  detailed  chemistry  model  (Sp)  and  the  overall 
reaction  model  (Sc)  .  The  percentage  error  in  the  one  reaction  model  is 
defined  as  UK)  ( S0  -  Sp)/Sp.  The  agreement  in  flame  speeds  is  excellent. 

Tables  1  and  4  show  the  corresponding  results  for  a  set  of  hydrogen/air 
flames.  Here,  the  overall  activation  energies  also  decrease  for  hotter  and 
faster  flames,  for  reasons  similar  to  those  discussed  for  the  ozone  case.  For 
a  slow  flame,  only  small  amounts  of  the  radical  species  (H,  OH,  0,  HO-))  are 
produced.  The  chain  branching  reactions  (H  +  "”011  +  0  and  0  +  If,  *  OH  + 

H),  which  have  relatively  high  activation  energies,  tend  to  be  rate 
limiting.  For  a  fast,  hot  flame,  large  amounts  of  radicals  are  created.  The 
flame  propagation  is  to  a  large  extent  determined  by  low  activation  energy 
reactions  (OH  +  1^  ►  Ho0  +  11). 

Results  for  methane/air  flames  are  given  in  Tables  r>  and  b.  As  before, 
the  overall  activation  energies  increase  for  slower  flames.  All  of  these 
activation  energies  are  relatively  high,  higher  than  for  any  of  the  reactions 
in  the  detailed  mechanism.  The  detailed  model  assumes  a  sequence  of  reactions 
( CH^  ►  CH^  *■  CH^O  *  CHO  *  CO  *  G^).  In  fitting  tin's  by  a  one  reaction 
model,  a  high  activation  energy  seems  to  best  fit  the  overall  rate  of  the 
sequence . 
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TABLF.  1.  ONK  REACTION  MODEL  PARAMETERS  FOR  OZONF./OXYGKN  FLAMES.* 


%  (>3 

CP 

Tb 

0 

pA 

A 

E 

20 

.244 

1094 

710 

7.32K-8** 

1 . 44  KH 

15.3 

2  3 

.248 

1256 

710 

7.20E-8 

5. 10E7 

12.8 

SO 

.259 

1943 

709 

6.86K-8 

1.39E7 

7.9 

75 

.265 

2425 

689 

6.71 K-8 

1 .33E7 

6.5 

100 

.269 

2702 

646 

6.63E-8 

1  .  59  F. 7 

6.0 

*  Units  are 

in  list  of 

symbols 

**  7.32K-8 

is  read  as 

7.32  x  10"8 

TABLE  2. 

A  COMPARISON  OF  FLAME 

SPEEDS  FOR  THE  DETAILED 

CHEMISTRY  MODEL  AND  THE  ONE 

REACTION  MODEL  FOR 

OZONE/OX YCEN  FLAMES. 

*  °3 

SD 

So 

Z  Error 

20 

33.1 

34.1 

3 

25 

61.5 

64.1 

4 

50 

241.6 

253.7 

5 

75 

391.1 

409.2 

5 

100 

496.8 

517.5 

4 

TABLE  3.  ONE 

REACTION 

MODEL  PARAMETERS  FOR  HYDROGEN/ A1 R  FLAMES 

tH-, 

CP 

Th 

0 

pA  A 

E 

15 

.297 

1473 

28890 

6.63E-8  1.58E17 

15. 

20 

.319 

1835 

28803 

6.60E-8  1.48E16 

9.i 

30 

.367 

2398 

27264 

6.72E-8  5.51E15 

7. 

40 

.416 

2214 

28604 

7.32E-8  2.06E15 

6 .' 

50 

.479 

1939 

28827 

7.65E-8  9. 61  El  4 

6 . 

60 

.567 

1642 

28890 

7.66K-8  4.84E14 

6. 

70 

.703 

1329 

28898 

7.32E-8  3.47E14 

8. 

75 

.803 

1166 

28890 

7.01E-8  7 .60814 

1  1  . 
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TABLE  4.  A  COMPARISON  OF  FLAME  SPEEDS  FOR  THE  DETAILED  CHEMISTRY 
MODEL  AND  THE  ONE  REACTION  MODEL  FOR  HYDROGEN/AIR  FLAMES. 


X  H2 

SD 

So 

%  ERROR 

15 

33.6 

32.9 

-2 

20 

100.5 

89.1 

-11 

20 

217.9 

203.8 

-7 

40 

280.6 

270.5 

-4 

50 

251.9 

251.9 

-0 

60 

178.0 

182.4 

3 

70 

82.4 

84.6 

3 

75 

35.5 

36.3 

2 

TABLE  5.  ONE  REACTION  MODEL  PARAMETERS  FOR  MKTHANK/AIR  FLAMES. 


':CH4 

cr 

Tb 

Q 

pX 

A 

E 

6.5 

.  300 

1778 

11960 

5.72E-8 

3. 37  E 1 8 

35.1 

7.5 

.308 

1960 

1 1908 

5 . 66E-8 

8. 15E17 

30.1 

8.5 

.316 

2121 

11768 

5.62E-8 

5.36E17 

28.1 

9.5 

.323 

2231 

11355 

5.61 E-8 

1 . 0 1  E  1 8 

29.1 

10.5 

.327 

2202 

11396 

5.68E-8 

9 . 1  5  E 1  7 

29.0 

1  1.5 

.330 

2107 

10002 

5.80E-8 

1.86E18 

31  .1 

12.5 

.334 

2009 

10582 

5.94E-8 

3 . 6  2  E 1 8 

33.8 

13.5 

.337 

1914 

10164 

6.I0E-8 

4. 10EI8 

35.6 

TABLE  6.  A  COMPARISON  OF  FLAME  SPEEDS  FOR  THE  DETAILED  CHEMISTRY 
MODEL  AND  THE  ONE  REACTION  MODEL  FOR  METHANE/ AIR  FLAMES 


%  ch4 

sn 

So 

Z  ERROR 

6.5 

16.0 

18.4 

15 

7.5 

25.3 

27.7 

9 

8.5 

33.8 

34.6 

2 

9.5 

39.8 

38.4 

-4 

10.5 

43.9 

40.8 

-7 

11.5 

42.6 

38.2 

-10 

12.5 

34.2 

30.1 

-12 

13.5 

23.5 

20.6 

-12 
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In  summary,  an  overall  reaction  model  appears  to  he  a  quite  accurate 
representat Ion  of  a  wide  variety  of  flames  and  stoichiometries.  The 
activation  energies  generated  are  relatively  low  compared  to  the  range  otten 
assumed  for  overall  reactions.  The  overall  reaction  normally  does  not 
correlate  with  any  single  reaction  in  the  detailed  model,  and  is  a  (unction  ot 
stoichiometry  as  well  as  the  detailed  kinetics. 

lip  to  this  point  we  have  arbitrarily  chosen  the  order  +  v^  (I  tor 
ozone,  3  for  hydrogen/air  and  methane/air),  based  on  the  assumed  overall 
reactions.  Some  results  were  computed  with  different  choices  lor  the  order 
but  with  the  ratio  Vp/v  constant.  The  consistent  result  was  that  as  the 
order  increases,  the  fitting  parameters  A  and  K  also  increased.  However,  the 
accuracy  of  the  fit  was  almost  independent  of  the  order  assumed. 

Also,  a  number  of  three  parameter  fits  were  made,  where  A,  and  K  were 

all  determined  by  a  least  squares  fit.  The  order  so  determined  varied 
substantially  with  stoichiometry.  The  accuracy  of  the  three  parameter  til  was 
only  a  slight  improvement  over  the  two  parameter  fit,  and  for  a  few  cases  it 
was  less  accurate. 

Thus,  over  a  large  interval,  the  order  can  he  arbitrarily  chosen  without 
affecting  the  accuracy  of  the  overall  reaction  fit.  An  increase  in  the  order 
also  increases  the  parameters  A  and  K.  Preliminary  work,  however,  indicates 
that  the  reaction  order  is  not  arbitrary  if  the  same  kinetics  are  applied  to 
different  pressures. 

VII r.  CONCLUSIONS 

The  aim  of  this  work  has  been  to  explore  the  applicability  of  the  overall 
reaction  rate  concept  to  premixed  laminar  flames  with  complex  chemistry.  The 
results  demonstrate  the  ability  of  a  single  reaction  formalism,  having  a  rate 
constant  of  simple  Arrhenius  form,  to  track  the  temperature,  heat  release,  and 
to  some  extent,  species  profiles  obtained  from  mul ti reaction  network 
calculations.  A  straightforward  method  using  a  least  squares  fitting 
procedure  is  presented  which  extracts  the  kinetic  parameters  directly  from  the 
heat  release  profile  as  a  function  of  temperature.  These  global  kinetic 
parameters  are  then  tested  for  their  ability  to  reproduce  the  overall  features 
of  the  complex  reaction  network,  thereby  allowing  a  quantitative  assessment  to 
be  made  of  their  utilization  as  a  contracted  representation  for  the 
thermochemistry  of  a  multicomponent  flame.  We  believe  our  approach  is  the 
first  attempt  to  derive  overall  reaction  rates  systematically  by  using  the 
"complete"  information  obtainable  from  a  detailed  laminar  flame  model.  This 
differs  significantly  from  previous  methods  which  were  tailored  to  extract 
kinetic  parameters  from  a  very  limited  amount  of  experimental  data. 
Consequently,  the  results  obtained  from  this  study  indicate  the  exact  extent 
to  which  our  global  reaction  model  is  capable  of  condensing  a  complex  reaction 
network  to  a  single  overall  reaction. 

One  can  organize  the  findings  into  several  echelons  of  detail.  The 
first,  and  most  detailed,  involves  how  well  a  single  formal  reaction  can  mimic 
the  net  rate  of  heat  release  (due  to  the  multireaction  network)  as  a  function 
of  position  in  the  flame.  A  related  issue  is  the  determination  of  the  overall 
reaction  parameters  which  best  represent  this  heat  release  profile.  At  an 
intermediate  level  of  detail  a  comparison  can  be  made  between  the  temperature 
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and  reactant  profiles  as  computed  by  the  multireaction  code  and  those 
predicted  by  the  overall  reaction.  Finally,  one  can  compare  the  flame  speed 
predicted  by  the  overall  reaction  model  with  the  multireact<on  code  value. 

The  levels  of  decreasing  detail  correspond  to  increasing  degrees  of 
integration  over  the  underlying  reaction  processes.  One  might  expect, 
therefore,  that  the  flame  speed  would  be  less  sensitive  than  the  temperature 
profile  to  imperfect  fidelity  in  the  overall  reaction  rate  profile.  This 
expectation  is  generally  borne  out  by  the  results.  Calculation  of  flame  speed 
with  the  overall  reaction  is  accurate  for  all  flames  to  about  }*>%  at  worst  and 
frequently  in  the  5%  range.  Accuracy  of  the  temperature  profiles  computed 
using  the  overall  reaction  is  generally  good  and  is  limited  primarily  by  our 
requiring  that  the  proper  adiabatic  temperature  be  achieved.  Attainment  of 
this  complete  heat  release  is  often  considerably  delayed  due  to  slow  radical 
recombination  in  the  late  stages  of  the  flame.  The  fidelity  with  which  the 
net  heat  release  rate  can  be  approximated  by  a  single  reaction  varies  somewhat 
with  the  flame  identity.  The  ozone  and  f^/air  representation  is  quite 
credible  over  a  very  wide  temperature  range,  whereas  for  the  CH^/air  the  range 
is  rather  narrow. 

Another  finding  of  this  study  is  that  the  kinetic  parameters  of  the 
overall  reaction  do  not  correlate  with  any  single  reaction  step  of  the  complex 
kinetics.  Furthermore,  the  values  of  these  parameters  change  with 
stoichiometry.  The  use  of  a  single  set  of  kinetic  parameters  to  represent  a 
number  of  different  stoichiometries  will  lead  to  further  compromises  in  the 
thermal  structure  description.  We  have  not  examined  this  trade-off 
exhaustively,  but  our  feeling  is  that  one  would  be  better  advised  to  treat  the 
overaLl  reaction  parameters  as  functions  of  stoichiometry. 

Preliminary  work  suggests  that  the  overall  reaction  concept  retains 
accuracy  over  a  wide  range  of  pressures  and  initial  temperatures.  The 
effective  values  of  the  kinetics  constants  vary  only  slightly  with  changes  in 
the  initial  temperature;  variation  with  pressure  is  slight  for  some  flames  and 
large  for  others.  These  considerations  will  be  addressed  in  a  subsequent 
paper . 

We  have  shown  that  an  overall  reaction  model  can  quite  accurately 
reproduce  many  of  the  features  of  a  detailed  chemistry  model  for  laminar 
flames,  despite  the  vast  simpl i f icat ion  involved.  The  extent  to  which  the 
overall  reaction  parameters  determined  from  this  problem  can  be  applied  to 
other  flow  conditions  remains  an  important  question  for  future  study. 
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LIST  OF  SYMBOLS 


A  =  Lumped  parameter,  cent imeter-gram-second  units. 


Ci 


-  concent  rat  ion  of  species  j,  Tiole-cm 


c*  * 


Dim 


specific  heat  of  the  mixture,  cal-g  ^-K 

=  diffussion  coefficient  of  species  i  in  the  mixture, 

2-1 

cm^-s 


K  = 


hi 


k  = 


mo  = 


activation  energy,  kcal-mole-*. 

-  specific  enthalpy  of  species  i,  eal-g  '. 
rate  coefficient  in  cent imeter-mole-second  units. 

e i  -2  -1 

mass  flux,  g-cm  -s  . 


M 


-1 


l 

P  = 
M  = 
Q  = 

R  = 

R 


=  molecular  weight  of  species  i,  g-mole  . 
pressure,  atm. 

_  O  _  I 

heat  release  per  unit  volume,  cal-cm  J-s  1  . 
heat  of  reaction  per  unit  mass  of  fuel,  c.al-g~* 
gas  constant  =  1.9872F.-3  kcal-mole~i-K~1  . 

=  gas  constant  =  82.05  cm^-atm-mole- . 


Ri  = 
S  = 


—  3  '1 

rate  of  production  of  species  i,  mole-cm  -s 
velocity  of  the  flame  relative  to  the  unburned  mixture, 


T  =  temperature,  K. 


,-l 


Yi 


fluid  velocity,  cm-s 
spatial  coordinate,  cm. 

=  mass  fraction  of  species  i. 

thermal  conductivity  of  the  mixture,  cal-cm-*-s-*-K-* . 
density  of  the  mixture,  g-cm-^. 
stoichiometric  coefficient. 
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Columbus,  OH  43201 

2  Exxon  Research  &  Engineering 
ATTN:  A.  Dean 

M.  Oon 
P.O.  Box  45 
Linden,  NJ  07036 

1  Ford  Aerospace  and 

Communications  Corp. 

D1VAD  Division 
Div.  Hq.,  Irvine 
ATTN:  D.  Williams 
Main  Street  6  Ford  Road 
Newport  Beach,  CA  92663 

1  General  Electric  Armament 

&  Electrical  Systems 
ATTN:  M.d .  Bui man 

Lakeside  Avenue 
Burlington,  VT  05402 

1  General  Electric  Company 

ATTN:  M.  Lap p 
Schenectady,  NY  12301 

1  General  Electric  Ordnance 

Systems 

ATTN:  .J.  Mandzy 
100  Plastics  Avenue 
Pittsfield,  MA  01203 

1  General  Motors  Rsch  Labs 

Physics  Department 
ATTN;  J.H.  Bechtel 
Warren,  Ml  48090 
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3  Hercules,  Inc. 

Alleghany  Ballistics  Lab. 
ATTN:  R.R.  Miller 
P.O.  Box  210 
Cumberland,  MD  21S01 

3  Hercules,  Inc. 

Bacchus  Works 
ATTN:  K.P.  McCarty 
P.0.  Box  98 
Magna,  UT  84044 

1  Hercules,  Inc. 

AFATL/DLDL 

ATTN:  R.L.  Simmons 

Eglin  AFB,  FL  32542 

[  Honeywell,  Inc. 

Defense  Systems  Division 
ATTN:  D.E.  Broden/ 

MS  MN50-2000 
bOO  2nd  Street  NE 
Hopkins,  MN  55343 

1  IBM  Corporation 

ATTN:  A.C.  Tam 
Research  Division 
5600  Cottle  Road 
San  Jose,  CA  95193 
1  Director 

Lawrence  Livermore  Nat ' 1  Lab 
P.0.  Box  808 
ATTN:  C.  Westbrook 
Livermore,  CA  94550 

1  Lockheed  Palo  Alto  Rsch  Lab 

ATTN:  George  Lo 
3251  Hanover  Street 
Dept.  52-35/B204/2 
Palo  Alto,  CA  94304 

1  Los  Alamos  National  Lab 

ATTN:  B.  Nichols 

T7.  MS-B284 
P.0.  Box  1663 
Los  Alamos,  NM  87545 


1  Science  Applications,  Inc. 

ATTN:  H.  S.  Pergament 
1100  State  Road,  Bldg.  N 
Princeton,  NJ  08540 

1  Olin  Corporation 

Smokeless  Powder  Operations 
ATTN:  R.L.  Cook 
P.O.  Box  222 
St.  Marks,  FL  32355 

1  Paul  Gough  Associates,  Inc. 

ATTN:  P.S.  Cough 

1048  South  Street 
Portsmouth,  NH  03801 

2  Princeton  Combustion 

Research  Laboratories 
ATTN:  M.  Summerfield 
N.A.  Messina 

475  US  Highway  One  North 
Monmouth  Junction,  NJ  08852 

1  Pulsepower  Systems,  Inc. 

ATTN:  L.C.  Elmore 
815  American  Street 
San  Carlos,  CA  94070 

1  Rockwell  International  Corp. 
Rocketdyne  Division 

ATTN:  J.E.  Flanagan/ HB02 
6633  Canoga  Avenue 
Canoga  Park,  CA  91304 

2  Sandia  National  Laboratories 
Combustion  Sciences  Dept. 
ATTN:  R.  Cattolica 

D.  Stephenson 
Livermore,  CA  94550 

1  Science  Applications,  Inc. 

ATTN:  R.B.  Edelman 
9760  Owensworth  Avenue 
Chatsworth,  CA  91311 
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1  Science  Applications,  Inc. 

ATTN:  R.  B.  Edelman 
9760  Owensworth  Avenue 
Chatsworth,  CA  91311 

1  Space  Sciences,  Inc. 

ATTN:  M.  Farber 
Monrovia,  CA  91016 

A  SRI  International 

ATTN:  S.  Barker 

D.  Crosley 
D.  Golden 
Tech  Lib 

333  Ravenswood  Avenue 
Menlo  Park,  CA  94025 

1  Stevens  Institute  of  Tech. 

Davidson  Laboratory 
ATTN:  R.  McAlevy,  III 

Hoboken,  NJ  07030 

1  Teledyne  Me . ormack-Selph 

ATTN:  C.  Leveritt 
3601  Union  Road 
Hollister,  CA  95023 

1  Thiokol  Corporation 

Elkton  Division 
ATTN:  W.N.  Brundige 
P.O.  Box  241 
Elkton,  MD  21921 

3  Thiokol  Corporation 

Huntsville  Division 
ATTN:  D.A.  Flanagan 
Huntsville,  AL  35807 

3  ThiokoJ  Corporation 

Wasatch  Division 
ATTN:  J.A.  Peterson 
P.O.  Box  524 

Brigham  City,  UT  84302 

l  United  Technologies 

ATTN:  A.C.  Eckbreth 
East  Hartford,  CT  06108 


2  United  Technologies  Corp. 
ATTN:  R.S.  Brown 

R.O.  McLaren 
P.O.  Box  358 
Sunnyvale,  CA  91086 

1  Universal  Propulsion  Company 

ATTN:  H..I.  McSpadden 
Black  Canyon  Stage  1 
Box  1140 

Phoenix,  AZ  85029 

1  Veritay  Technology,  Inc. 

ATTN:  E.B.  Fisher 
P.O.  Box  22 

Bowmansville,  NY  14026 

1  Brigham  Young  University 

Dept,  of  Chemical  Engineering 
ATTN:  M.W.  Beckstead 
Provo,  UT  84601 

1  California  Institute  of  Tech. 

Jet  Propulsion  Laboratory 
4800  Oak  Grove  Drive 
Pasadena,  CA  91103 


1  California  Institute  of 

Technology 

ATTN:  F.E.C.  Culick/ 

MC  301-16 
204  Karman  Lab. 

Pasadena,  CA  91125 

1  University  of  California, 

Berkeley 

Mechanical  Engineering  Dept. 
ATTN:  J.  Daily 
Berkeley,  CA  94720 

1  University  of  California 

Los  Alamos  National  Lab. 

ATTN:  T.D.  Butler 

P.O.  Box  1663,  Mail  Stop  B2 1 6 

Los  Alamos,  NM  8764 S 
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2  University  of  California, 

Santa  Barbara 
Quantum  Institute 
ATTN:  K.  Schoftel.l 

M.  Steinberg 

Santa  Barbara,  CA  S310b 

1  University  of  Southern 

California 
Dept,  of  Chemistry 
ATTN:  S.  Benson 
Los  Angeles,  CA  90007 

l  Case  Western  Reserve  Univ. 

Div.  of  Aerospace  Sciences 
ATTN:  .1.  Tien 
Cleveland,  OH  44135 


]  Hughes  Aircraft  Company 

ATTN:  T.E.  Ward 
8433  Fall  brook  Avenue 
Canoga  Park,  CA  9130.5 

1  University  of  Illinois 

Department  of'  Mech  Png 
ATTN:  H.  Krier 

M1-B,  120b  IV  Croon  Street 
Urbana,  IL  81801 

1  Johns  Hopkins  Univers i  tv/ API. 

Chemical  Propulsion 
Information  Agency 
ATTN:  T.W.  Christian 
Johns  Hopkins  Road 
Laurel,  MD  207(17 


1  Cornell  University  1  University  of  Minnesota 

Department  of  Chemistry  Dept,  of  Mechanical 

ATTN:  E.  Grant  Engineering 

Baker  Laboratory  ATTN:  K.  Fletcher 

Ithaca,  NY  14853  Minneapolis,  MN  55455 


1  Univ.  of  Dayton  Rsch  Inst. 

ATTN:  D.  Campbell 

AFRPL/PAP  Stop  24 
Edwards  AFB,  CA  93523 

1  University  of  Florida 
Dept,  of  Chemistry 
ATTN:  J.  Winefordner 
Gainesville,  FL  32bll 

3  Georgia  Institute  of 

Technology 
School  of  Aerospace 
Engineering 
ATTN:  E.  Price 
Atlanta,  GA  30332 

2  Georgia  Institute  of 

Technology 
School  of  Aerospace 
Engineering 
ATTN:  W.C.  Strahle 
B.T.  Zinn 

Atlanta,  GA  30332 


4  Pennsylvania  State  University 

Applied  Research  Laboratory 
ATTN:  G.M.  Faeth 
K.K.  Kuo 

H.  Palmer 
M.  Micci 

University  Park,  PA  lbS02 

1  Polytechnic  Institute  of  NY 

ATTN:  S.  Lederman 

Route  110 

Farmingdale,  NY  11735 

2  Princeton  University 

Forrestal  Campus  Library 
ATTN:  K.  Brezinsky 

I.  Classman 
P.0.  Box  710 
i.in.eton,  NJ  08540 

1  pri>  eton  University 

MAE  r'ep.  . 

ATTN:  F.A.  Williams 

Princeton,  NJ  08511 
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Purdue  University 

1 

Stanford  University 

School  ot  Aeronautics 

Dept,  of  Mechanical 

and  Astronautics 

Eng i nee r i ng 

ATTN:  R.  Click 

ATTN:  R.  Hanson 

J.R.  Osborn 

Stanford,  CA  91.308 

Crissom  Hall 

West  Lafayette,  IN  47907 

2 

University  of  Texas 

Purdue  University 

Dept,  of  Chemistry 

ATTN:  W.  Cardiner 

School  ot  Mechanical 

H.  Schaefer 

Eng i in  ering 

Austin,  TX  78712 

ATTN:  N.M.  I,aurendeau 

S.N.B.  Murthy 

1 

University  of  Utah 

D.  Sweeney 

Dept,  of  Chemical  Engineering 

TSPC  Chaffee  Hall 

ATTN:  0.  Flandro 

West  Lafayette,  IN  47906 

Salt  Lake  City,  1'T  84112 

Rensselaer  Polytechnic  Inst. 

1 

Virginia  Polytechnic- 

Dept,  of  Chemical  Engineering 

Institute  and 

Ai'TN:  A.  Fontijn 

State  University 

Troy,  NY  12181 

ATTN:  d.A.  Sehetz 

Southwest  Research  Institute 
ATTN:  Robert  F..  White 

Aberdeen 

Blacksburg,  VA  24061 

Proving  Ground 

A.B.  Wenzel 

San  Antonio,  TX  ”8228 

Dir,  USAMSAA 

ATTN:  DRXSY-I) 

DKXSY-MF ,  H.  Cohen 
Cdr ,  USATK.COM 

ATTN:  DRSTK-TO- F 

Dir,  USACSL,  Bldg.  K331b,  KA 
ATTN:  DRDAR-CI.B-PA 

DRDAR-CLN 
DRDAR-OU-L 


uslr  l valuation  or  rlport 


Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  place 
in  the  mail.  Your  comments  will  provide  us  with  information  for 
improving  future  reports. 

1 .  BRL  Report  Number _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.)  _  _  _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) _ 


b.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name : _ 

Telephone  Number; _ 

Organization  Address:  _ _ 
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